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CHAPTER  1 


INTRODUCTION 

The  problem  of  electromagnetic  coupling  from  one 
region  to  another  region  through  an  aperture  in  a 
conducting  plane  can  be  formulated  in  terms  of  two 
admittance  operators,  one  for  each  region  Cl  3.  These 
admittance  operators  are  complex  and  symmetric.  Recently, 
a  theory  for  the  characteristic  modes  for  apertures  has 
been  proposed  C23.  The  characteri Stic  currents  are 
defined  as  the  eigenfunctions  of  a  certain  generalized 
eigenvalue  equation  involving  the  admittance  operators  and 
their  real  parts.  Because  of  the  particular  choice  of  the 
eigenvalue  equation,  the  characteristic  currents  are  real 
(or  equiphase)  and  orthogonal  with  respect  to  the 
admittance  operator,  its  real  part,  and  its  imaginary  part 
over  the  aperture.  Furthermore,  the  characteri sti c  fields 
produced  by  the  characteristic  currents  are  orthogonal 
over  the  radiation  sphere.  For  small  apertures,  the 
characteri stic  mode  theory  reduces  to  an  augmented  Bethe 
hole  theory,  i.e.,  the  aperture  is  described  by  a 
susceptance  term  related  to  the  polarizability,  plus  a 
conductance  term. 

Of  aperture  problems,  the  problem  of  an  infinitely 
long  slot  in  a  conducting  plane  (see  Figure  1-1)  has  been 


igure  1-1.  An  infinitely  long  slot  in  a 
conducting  plane. 


considered  most  over  the  years  since  it  was  first 
addressed  by  Lord  Rayleigh  L33 .  In  his  solution.  Lord 
Rayleigh  relied  on  reducing  the  problem  to  known 
solutions  of  hydrodynamic  or  electrostatic  problems,  but 
treated  only  the  special  case  of  a  sufficiently  narrow 
slot  (a  slit)  C4,  Page  2733.  The  narrow  slot  has  since 
continued  to  receive  particular  attention.  Morse  and 
Rubenstein  C53  provided  an  exact  solution  for  the  slit. 
Their  analysis  was  based  on  solving  the  Helmholtz  equation 
in  elliptic  coordinates  using  separation  of  variables, 
and,  as  a  result,  involved  series  of  Mathieu  functions. 
Later,  Barakat  C63  generalized  their  work  to  a  slot  in  a 
conducting  plane  separating  mediums  with  different 
electromagnetic  properties.  Sommer f eld  [4,  Section  393 
reduced  the  associated  boundary  value  problem  into  an 
integral  equation  for  the  electric  field  in  the  slot  for 
uniform  transverse  electric  and  magnetic  plane  waves 
normally  incident  on  the  slot.  Although  the  formulation 
is  general  and  valid  for  wide  slots  of  arbitrary  width, 
the  solution  was  given  only  under  narrow  slot 
approximations.  Millar  E73,  extending  Sommerfeld’s 
equations  to  obliquely  incident  plane  waves,  solved  for 
the  electric  field  in  the  slot  in  the  form  of  a  power 
series  in  the  ratio  of  slot  width  to  wavelength.  Houlberg 
CS3  then  generalized  the  work  of  Millar  to  the  two-medium 


problem. 


Variational  techniques  C9,  Section  7-121  also 


■found  their  application  to  the  solution  of  narrow  slots. 
More  recently,  Butler  and  Milton  CIO!  used  a  series 
expansion  of  Chebyshev  polynomials  to  obtain  exact 
solutions  for  the  equivalent  magnetic  current  on  the 
narrow  slot  for  a  general  excitation. 

Effective  solutions  for  wide  slots  of  width 
comparable  to  the  wavelength,  however,  have  only  become 
possible  with  the  advance  of  solution  techniques  based  on 
the  method  of  moments  till.  Typically,  an  operator, 
usually  integral,  equation  is  derived  for  some  unknown 
parameter  from  which  all  other  quantities  can  be  derived. 
As  in  narrow  slot  solutions,  the  unknown  parameter  is 
approximated  by  a  linear  combination  of  some  known 
functions,  but,  in  contrast  to  narrow  slot  solutions,  the 
kernel  of  the  integral  equation  is  not  replaced  by  its 
small  argument  approximation.  The  inner  products  of  the 
integral  equation  with  testing  functions  form  a  matrix 
equation  whose  solution  determines  the  unknown 
coefficients.  Butler  and  Umashankar  C123  utilized  the 
generalized  admittance  formulation  to  solve  the 
two-medium  problem.  Chou  and  Adams  C133  treated 
single  and  double  slots  in  an  unbounded  medium  by 
converting  the  slot  problem  to  that  of  a  strip  using 
Babinet's  principle  C9,  Section  7-123.  Other  varieties  of 
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slot  problems  were  also  considered.  Lewis  E143  studied 
the  effect  of  covering  the  slot  by  a  dielectric  sheath, 
while  Nevels  and  Butler  E153  solved  the  problem  o-f  a  slot 
in  a  ground  screen  covered  by  a  dielectric  slab. 
Electromagnetic  transmission  through  a  filled  slot  in  a 
thick  conducting  plane  was  investigated  by  Auckland  and 
Harrington  C163,  L171. 

In  this  report,  the  theory  of  characteristic  modes 
for  apertures  is  specialized  to  an  infinitely  long  slot  in 
a  conducting  plane  in  an  unbounded  medium.  The  theory  is 
then  applied  to  a  complete  solution  of  the  slot  problem 
when  the  slot  is  illuminated  by 

1-  a  uniform  transverse  electic  (to  the  slot  axis/ 
plane  wave, 

2-  a  uniform  tranverse  magnetic  (to  the  slot  axis) 
plane  wave. 

Specifically,  the  characteristic  currents  and  field  modes 
are  computed  for  different  slots.  These  are  then  used  to 
compute  the  quantities  and  parameters  of  importance 
usually  encountered  in  electromagnetic  field  compatibility 
problems,  such  as  the  equivalent  magnetic  current  of  the 
slot,  the  transmission  coefficient,  and  transmitted  field 
pattern  far  from  the  slot.  Analytic  expressions  for  the 
special  case  of  the  narrow  slot  are  also  given.  It  is 


m 

shown  that  the  theory  of  character! stic  modes  is  a  very 


useful  theoretical  tool.  It  allows  -for  a  modal  expansion 
o-f  all  quantities  o-f  interest,  and  can  be  applied  to 
narrow  slots  as  well  as  to  general  slots  o-f  arbitrary 
width.  Furthermore,  it  is  shown  that  it  is 
computationally  very  efficient.  It  exhibits  -features  and 
speed  o-f  convergence  that  make  it  more  attractive  than 
many  other  numerical  solutions. 

The  theory  of  character i sti c  modes  is  applied  in  the 
report  only  to  a  slot  in  a  conducting  plane  in  an 
unbounded  medium.  Such  a  restriction  is  not  essential, 
and  is  intended  only  to  simplify  the  theory.  The 
extension  to  the  two-medium  problem  is  straightf orward, 
and  can  be  accomplished  with  a  few  obvious  changes  in  the 
present  analysis.  The  conclusions  drawn  here  for  the  one- 
medium  problem  are  expected  to  hold  for  the  two— medium 
problem.  Finally,  a  word  about  the  organization  of  the 
report  is  in  order.  In  Chapter  2,  the  theory  of 
characteri sti c  modes  for  a  slot  in  a  conducting  plane  in 
an  unbounded  medium  is  presented.  In  Chapters  3  and  4, 
the  theory  of  characteri sti c  modes  is  applied  to 
transverse  electric  and  magnetic  plane  wave  excitations  of 
the  slot,  respectively.  In  Chapter  5,  conclusions  are 
drawn  and  some  final  remarks  are  given. 


V, 


CHAPTER  2 


THE  CHARACTERISTIC  MODES  FOR  A  SLOT 
IN  A  CONDUCTING  PLANE 

The  character! sti c  currents  and  -fields  for  a  slot  in  a 
conducting  plane  in  an  unbounded  medium  are  defined  in  this 
chapter.  The  characteristic  currents  are  the  solutions  of 
an  eigenvalue  equation  relating  the  tangential  components  of 
the  magnetic  field  at  the  slot.  The  characteristic  fields 
are  the  fields  produced  by  the  characteristic  currents  on 
the  slot  in  the  presence  of  the  complete  screen.  These 
currents  and  field  modes  are  then  used  for  the  solution  of 
the  slot  problem.  A  Gal er kin  solution  of  the  eigenvalue 
equation  is  given  at  the  end  of  the  chapter. 

2- 1  Derivation  of  the  Operator  Equation 

Let  the  excitation  of  the  slot  be  a  field  (E1 ,  H/ ) 
incident  from  the  left  of  the  screen  (z<0) .  This  field  is 
called  the  incident  field.  It  is  the  field  that  would  exist 
if  the  conducting  screen  were  absent.  Because  of  the 
presence  of  the  screen  with  the  slot,  part  of  the  incident 
field  is  reflected  back,  while  the  rest  is  transmitted  into 
the  z >0  half-space.  The  total  field,  incident  plus 
scattered,  <E,  H) ,  must  have  zero  electric  field  component 
tangent  to  the  screen,  and  continuous  tangential  electric 

Bel ow. 


and  magnetic  fields  across  the  slot. 


a  field 


equivalence  theorem  is  used  to  divide  the  problem  into  two 
decoupled  parts. 

Let  the  exciting  -field  be  incident  while  the  slot  is 
covered  by  a  perfect  conductor.  This  field,  often  referred 
to  as  the  short-circui t  field,  is  denoted  (E«=,  H-*=).  By 
the  field  equivalence  theorem  C9,  Section  3-5],  the  field  to 
the  left  of  the  screen  is  identical  with  (E**,  H“c)  plus  the 
field  (E(M) ,  H(M> >  produced  by  the  magnetic  current  sheet 

M  *  u»  x  E|  (2-1) 

>z=0 

on  the  slot  while  it  is  covered  by  a  perfect  conductor.  The 
field  to  the  right  of  the  screen  is  then  identical  with  the 
field  (E(-M),  H(-M) )  produced  by  the  magnetic  current  sheet 

-M  on  the  slot  while  it  is  covered  by  a  perfect  conductor. 
Figure  2-1  shows  the  equivalent  situations. 

The  total  tangential  electric  field  clearly  vanishes  at 
the  conducting  screen.  Furthermore,  the  tangential  electric 
field  is  continuous  across  the  slot  by  virtue  of  placing 
magnetic  current  sheets  of  opposite  signs  on  the  opposite 
sides  of  the  slot.  The  continuity  of  the  magnetic  field 
across  the  slot,  however,  requires  that 

( M )  +  H*-«=  =  H„(-M>  =  -  H*(M>.  (2-2) 

In  (2-2),  the  subscript  t  refers  to  the  tangential 
component,  and  the  last  equality  follows  because  of  the 


Original  problem  Equivalent  situation  Equivalent  situation 

for  e<0  for  *>0 

Figure  2-1.  The  equivalent  situations. 


linearity  of  the  magnetic  field, 
becomes 


Rearranging,  (2-2)  then 


& 


I 


*5 


-2  H*(M)  =  H*-=. 


(2-3) 


In  the  next  section,  the  characteristic  currents  of  the 
slot  are  defined.  These  currents  are  then  utilized  to  solve 
(2-3)  for  M. 


2-2  The  Characteristic  Currents  of  the  Slot 


The  operator  equation  (2-3)  can  be  rewritten  in  the 


form 


Y (M)  =  I 


(2-4) 


where 


Y(M)  =  -2  Ht (M) 


(2-5) 


I  =  . 


Since  M  and  I  have,  respectively,  the  dimensions  Volt/meter 
(V/m)  and  Amp&re/meter  (A/m),  Y  is  an  operator  with  the 
dimension  of  an  admittance. 


*  * 

Define  6  =  (Y  +  Y  ) /2  and  B  =  (Y  -  Y  )/(2j),  then 


Y(M>  =  6(M>  +  jB(M) 


(2-6) 


In  (2-6),  G  is  a  conductance  operator  and  B  is  a  susceptance 
operator.  Following  Harrington  and  Mautz  C2J,  the 


characteristic  currents  of  the  slot  are  defined  to  be  the 


eigenfunctions  Mr,  of  the  eigenvalue  equation 

Y  (Mr.)  =  y„  G  ( Mr. )  (2-7) 

so  normalized  that 

<Mr.,  G (Mr, )  >  =1.  (2-8) 

In  (2-8),  <-,  ->  denotes  the  inner  product 

+w 

<C,  D>  =  I  C  (x>  .  D(x)  dx  (2-9) 

— w 

* 

where  C  (x)  is  the  complex  conjugate  of  C(x)  and  ±w  are  the 
x  coordinates  of  the  edges  of  the  slot.  Put 

y„  =  1  +  jbn.  (2-10) 

Then,  using  (2-6)  and  (2-10),  (2-7)  becomes 

B(Mr»)  *  br.  G  (Mr.  ) .  (2-11) 

The  operators  G  and  B  are  self-adjoint,  since  for  any 
M»  and  M* 

* 

,  G(Ma>>  =  <Mi,  *(Y  +  Y  )  (M»)  > 

* 

■  <!*(Y  +  Y  >  (Mi),  Ma> 


=  <G(M») ,  M»> 


(2-12) 


til 


& 
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and  similarly  for  B.  Furthermore,  S  is  positive  definite, 
since  the  time-average  power  (per  unit  length  in  the  y 
direction)  radiated  into  the  z  >0  half-space 


P r  «d  —  *4  Re 


+w 

(I 


E  x  H  I  .  u*  dx 

U=0 


-w 


+w 


=  >i  Re 


£  M  .  H  dx  ] 


-w 


=  '4.  Re£  \  <  Y  (M)  ,  M>  j 


=  %  <G(M) ,  M> 


(2-13) 


is  always  positive.  The  second  and  the  last  equalities  in 
(2-13)  follow  from  (2-1)  and  (2-12)  respecti vel y.  It  is 
then  a  standard  practice  C18,  Section  1-253  to  prove  that 
all  b„,  and  hence  M„,  are  real,  and  that  can  be  chosen  to 
satisfy  the  orthogonal i ty  relationships 

1.  G(M„)>  = 


2.  <Mm,  B(M„)>  =  b„  Sn« 


3.  <Mm,  Y  (Mr*)  >  *  (1  +  jb„)  S, 


(2-14) 


where  is  the  Kronecker  delta  function  (0  if  m£n,  and  1 

i f  m=n ) . 


A  -V  * 


All  the  currents  on  a  slot  in  a  conducting  plane  in  an 
unbounded  medium  are  required  to  radiate  some  power  however 
small.  As  can  be  seen  -From  (2-11)  and  (2—13),  the 
characteristic  currents  corresponding  to  very  large  bn  are 
basically  non-radiating.  It  is  shown  in  later  chapters 
that,  when  the  slot  is  very  narrow,  all  the  characteristic 
currents  have  very  large  eigenvalues  Cb„>  and,  there-fore, 
are  essentially  non-radiating.  In  any  case,  all  currents 
are  required  to  exhibit  the  edge  property  (M„  =  CM-JW2  -  V*) 
and  M*,  =  0(1/ -Jw*  —  x * ) ,  as  x-f±w)  C19,  Section  1-43. 

2-3  Modal  Solution  of  the  Operator  Equation 

A  modal  solution  of  (2-4)  for  the  magnetic  current  M 
over  the  slot  is  obtained  in  this  section.  Put 

M  =  I  V„  M„  (2-15) 

n 

where  Ml,  are  the  characteristic  currents  of  the  slot,  and 
V„  are  complex  coefficients  to  be  determined.  Substituting 
(2-15)  into  (2-4),  it  then  becomes 

E  V„  Y (Mr,)  =  I.  (2-16) 

n 

Taking  the  inner  product  of  (2-16)  with  each  M*,,  there 


results 


14 


£  v„  <!im,  Y(M„)>  =  <Mm,  I>. 
n 


(2-17) 


$ 

M 


Because  of  the  third  orthogonal i ty  relationship  of  (2-14), 
only  the  mth  term  in  the  summation  survives.  Thus 


<M„,  Y(M*)>  =  <Mm,  I>. 


(2-18) 


Hence 


M 


vm  = 


<M«,  I> 


1  +  Jb, 


(2-19) 


Substituting  (2-19)  into  (2-15),  it  becomes 


<M«,  I> 

M  =  E  -  M«. 

n  1  +  jb„ 


(2-20) 


The  magnetic  current  M  given  by  (2-20)  is  called  the 
modal  solution  of  (2-4). 


2-4  Power  Consideration 


The  total  complex  power  entering  the  slot  is  basically 


Pin  =  ^  I  E  x  H  I  .  u.  dx. 

12=0- 

-w 


(2-21 ) 


Using  (2-1),  (2-4),  and  (2-5),  (2-21)  becomes 


=  'i  <1,  M> 

=  *  <Y(M)  ,  M>.  (2-22) 

Furthermore,  using  (2-20) ,  (2-22)  becomes 

- 

I<I,  M„ >  I a  .  (2-23) 


1  -  jb„ 

P,n  =  %  I  - 

n  1  +  b„a 


Since  the  time-average  power  (2-13)  radiated  into  the  z >0 
halt -space  is  equal  to  the  time-average  power  entering  the 
slot,  ?«B(M),  M>  is  equal  to  the  real  part  o-f  the  right-hand 
side  o-f  (2-23). 

A  parameter  sometimes  used  to  express  the  transmission 
characteristics  of  the  slot  is  the  transmission  coefficient 
T.  By  definition,  the  transmission  coefficient  of  the  slot 
is  the  ratio  of  the  time-average  power  transmitted  through 
the  slot  to  that  incident  on  the  slot  C9,  Section  7-123. 
Using  (2-23),  T  is  readily  found  as 

1  1 

T  =  —  Z  -  I<I,  h„>l2  (2-24) 

4P  n  1  +  b„* 

where  P  is  the  time-average  power  incident  on  the  slot. 

2-5  The  Characteri stic  Fields  of  the  Slot 

The  fields  (E„,  H„)  produced  by  the  characteri sti c 
currents  Mr,  are  called  the  characteri  sti  c  fields  of  the 


slot.  Orthogonality  relationships  -for  the  character! stxc 
fields  over  the  radiation  cylinder  can  be  obtained  -from 
those  for  the  characteristic  currents  by  means  of  the 


complex  Poynting  theorem  C9,  Section  1-103.  These 

relationships  are  dual  to  those  for  the  characteristic 
fields  of  the  conducting  body,  and  can  be  derived  in  a 
similar  manner  C203.  Thus 


1.  — 

n 


Em  .  En  dr  Sr 


Hm  .  H„  dr  =  Sr 


(2-25) 


Cy» 


3.  w 


* 

( 


Hn  -c  E»  .  E„)  dv  =  -b„  Sr 


in  duality  with  the  orthogonality  relationships  for  the 
conducting  body.  In  (2-25),  n  is  the  wave  impedance  of  the 
medium,  is  the  radiation  cylinder,  and  the  integration 
in  (2-25.3)  is  over  the  whole  space.  Figure  2-2  shows  the 
integration  domain  for  (2-25) . 

The  third  orthogonality  relationship  of  (2-25)  states 
that  the  difference  between  the  magnetic  and  electric  energy 
stored  in  any  characteristic  field  is  -b„^2co)  Joules  for 


every  one  Watt  of  radiated  power.  The  character i sti c  fields 


corresponding  to  positive  (negative)  b  have  predominantly 


stored  electric  (magnetic)  energy,  and  are  therefore 
referred  to  as  the  capacitive  (inductive)  modes  of  the  slot. 
The  characteri stic  fields  corresponding  to  b=0  are  then 
c a 1 1 ed  the  resonant  modes  of  the  slot . 

A  modal  expansion  of  the  field  radiated  by  M  in  terms 
of  the  characteristic  fields  £„  and  H„  can  be  readily 
obtained  using  the  modal  expansion  (2-15)  of  M: 


E(M)  =  £  V„  E (Mr, )  =  £  V„  E„ 
n  n 

(2-26) 

H(M)  =  £  V„  H (Mr, )  =  £  Vr,  Kr» . 
n  n 

In  Appendix  A,  it  is  shown  that  the  field  (E(M),  H(M) )  given 
by  (2-26)  converges  in  a  least-squares  sense  on  the 
radiation  cylinder. 


2-6  Solution  of  the  Eigenvalue  Equation 

An  exact  solution  of  the  eigenvalue  equation  (2-11)  for 
the  character i stic  currents  is  rather  difficult,  if  at  all 
possible.  An  approximate  solution  has  then  to  be  sought. 


Put 


9 


N 

S  =  LJ  S*  (2-27) 

k=  1 

L 

=  Z  U„K  f  „.  (2-28) 

k=l 

In  (2-27),  S*  are  non-overlapping  intervals  such  that 

N 

2w  =  Z  ( S*  I  (2-29) 

k=l 

where  IS*  I  is  the  length  of  the  kth  interval  (see  Figure  2- 
3).  In  (2—28),  each  f*«  is  a  real  -function  that  is  de-fined 
on  a  part  o-f  the  partition  and  vanishes  on  the  remainder  o-f 
it,  and  Unk  are  real  coefficients  to  be  determined. 
Substituting  (2-28)  into  (2-11),  it  becomes 

L  L 

I  §<!*>  =  b„  Z  IU*  G(f id  +  R  (2-30) 

k=l  k=l 

where  R  is  a  residual  term. 

A  Galerkin  solution  Cll,  Section  1-33  of  (2-11)  can  be 
obtained  by  requiring  that  R  be  orthogonal  to  all  f*,  viz.  » 

<fi ,  R>  =  0  1=1,2, ...,L.  (2-31) 

Thus,  taking  the  inner  product  of  (2-30)  with  each  f» ,  and 
enforcing  the  Galerkin  condition  (2-31),  there  then  results 


L 

I  U„*  <ix, 
k=l 


b 


L 

E  Unt< 
k=l 


G  ( -f  k  )  > 


(2-32) 


B(f  *)  >  = 


<f 


.1  > 


1  =  1,2 . L. 


In  matrix  -form,  (2-32)  becomes 


=  -»  =  4 

B  U„  =  bn  G  U„ 


(2-33) 


where  G  and  B  are  the  L  by  L  matrices 


G  =  CGx„]  =  C<fx,  G(1k)>D 

B  =  IBxul  =  C<f1?  B ( £ k )  >1 

and  U„  is  the  L  by  1  vector 

U„  =  L  U„  k  3  . 

The  constraint  equation  (2-8)  now  becomes 

->  T  =  -» 

U„  G  U„  =  1 

where  the  superscript  T  denotes  vector  transpose. 


(2-34) 


(2-35) 


(2-36) 


The  solution  of  (2-33)  determines  in  a  Galerkin  sense 
the  first  L  character i sti c  currents  of  the  slot. 


i&k 
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CHAPTER  3 


TRANSVERSE  ELECTRIC  (TE)  CASE 

In  this  chapter,  the  characteristic  currents  and 
magnetic  ‘ ields,  the  equivalent  magnetic  current,  and  the 
radiation  pattern  are  computed  for  a  general  slot  of 
arbitrary  width  in  a  conducting  plane  when  it  is 
illuminated  by  a  uniform  transverse  electric  (to  the  slot 
axis)  plane  wave.  Analytic  expressions  for  the  narrow  slot 
are  also  given. 


3-1  Basic  Formulation 

Let  a  plane  wave  be  incident  on  the  slot  at  an  angle  © 
from  the  left  of  the  screen  (see  Figure  3—1).  This  wave  is 
assumed  uniform  and  transverse  electric  to  the  y— axis,  and 
therefore  has  the  field  distribution 


-jk(x  sin©  +  z  cos©) 

E1  =  n Ceos©  u„  -  sine  u«3  e 

(3-1) 

-jk(x  sin©  +  z  cos©) 

H*  =  e  uv . 


Since  the  slot  is  uniform  along  the  y-axis,  and  since 
the  incident  magnetic  field  has  only  an  Hy.  component  that 
does  not  vary  with  y,  so  does  the  scattered  magnetic  field. 


2 


I 

i 

( 

i 

I 

Consequently,  the  total  -field  is  independent  o-f  y  and 
transverse  electric  to  the  y-axis.  It  then  -follows  from 
(2-1)  that  M  has  only  a  y-component  that  does  not  vary  with 

y: 


M  =  M  (x )  uv.  (3-2) 

The  field  to  the  left  of  the  screen  is  thus  given  by 
C9,  Section  3-123 


E(M>  =  -  V  x  F(M) 


k 

H  (li)  =  -j-  F  (M) 

-  n  - 


(3-3) 


plus  the  short-circuit  field 


=  -2nC jcose  sin(kz  cosd)  u^  +• 


-jkx  sin© 


sin©  cos(kz  cos©)  uf]  e 


(3-4) 


-jkx  sin© 

H““  =  2  cos(kz  cosd)  e  uv. 

In  (3-3),  F(M)  is  the  electric  vector  potential  produced  by 
M  xn  the  presence  of  the  complete  screen,  viz. 3  C9,  Section 
5-7  3: 


+w 


F  (M) 


1  r  (2) 

=  —  M  (x  ’  )  Ho  (kJz  = 

2j  J 

-w 


+  (x-x '  )*)  dx  •’  Uy 


(3-5) 


25 


(2) 

where  H©  is  the  Hankel  function  of  the  second  kind  and 
zero  order.  The  field  to  the  right  of  the  screen  is  then 
the  negative  of  that  in  (3-3). 

Substituting  (3-3),  (3-4),  and  (3-5)  into  (2-5),  Y ( M) 

and  I  are  readily  found  as 


+w 


Y  (M)  =  - 


(2) 

M(x')  H0  (klx-x’l)  dx’  u. 


-w 


-jkx  sin© 

I  =  2  e  uy. 


Since  C9,  Appendix  D3 


(2) 

Ho  (s)  =  J0(s)  -  j  N0  <  s) 


(3-6) 


(3-7) 


where  J0  and  No  are  the  Bessel  functions  of  the  first  kind 
and  zero  order,  and  of  the  second  kind  and  zero  order, 
respectively,  the  operators  G(M)  and  B(M)  in  (2-6)  are 
given  by 


.W, 

iUW 


Wf 

m 

wM 

iiwl 

m 

iVtl 


iuyuu 


G(M)  =  -  M  (x  )  Jo  <k I x-x ’  I )  dx ’ 


(3-8) 


B(M) 


■ 

k  r 

= -  M(x’)  N< 

n  J 


(k  ( x-x  ’ I >  dx 


In  the  next  section,  the  eigenvalue  equation  (2—11)  is 


solved  for  the  special  case  of  the  narrow  slot.  Thi 


numerical  solution  o-f  (2-11)  -for  the  general  slot  is  then 


consi dered. 


3-2  The  Narrow  Slot 


An  important  case  that  requires  special  consideration 


is  that  o-f  the  narrow  slot  (2kw<<l). 


Since 


k I x-x  *  •  <  2kw  <<  1 


(3-9) 


the  Bessel  -functions  J0  and  No  can  be  replaced  by  their 


small  argument  approximations  C9,  Appendix  D3: 


W 


1 


:!S 


a 


$ 

k 

$ 

Pa! 

c 


‘f 

*'l 

,r» 

t 

s 

i1 

V 


Jo  (k  I x-x '  I  )  ~  1 - (x  — x '  ) 

4 


2  rYk 

No  (klx-x’l)  ~  -  log  —  I  x— x  ’  I  j 


(3-10) 


where  Mlog“  denotes  the  natural  logarithm  and  Y  =1 . 7810724. 
6(H)  and  B(M)  are  then  given  by 


k 

r 

k  1 

i  G(M)  ~  — 

M(x’ )  : 

L - (x-x’)=* 

n  ‘J 

L 

4  -i 

(3-11) 


B(M) 


2k  P  fYk  , 

»  -  —  j  H(x')  log^— fx-x' I  j  dx' 


Furthermore,  I  does  not  vary  appreciably  over  S,  and  can 
therefore  be  approximated  by  the  first  two  terms  of  its 
Taylor  expansion  about  x=0.  That  is, 


I  s  2  -  2jkx  sinB  x€S 

or,  on  using  (3-4), 

k 

I  ~  Hv— (0)  +  j-  E«-=<0)  x  x<ES 

n 

where  (0)  is  written  for  the  point  (x,z)=(0,0>. 


(3-12) 


(3-13) 


The  characteristic  currents  of  the  narrow  slot  can  be 
determined  by  solving  (2-33)  with  S(M)  and  B(M)  given  by 
(3-11).  Since  the  slot  is  narrow,  however,  partitioning  o-f 

=  S. 

it  is  not  necessary.  The  lkth  elements  o-f  6  and  B  then 
become 


(x)  dx  fi«(x 


’)  j^l - (x-x’)*j  dx’ 


(3-14) 


2k  I*  f  fYk 

=  -  —  J  fi  (x)  dx  J  f*(x’  >  logj^ — Ix-x’IJ  dx '  . 


An  appropriate  choice  o-f  fi«  that  satisfies  the  edge 
requirements  on  Mv,  and  is  compatible  with  (3-13),  is 


f*(x)  = 


k=l . 2. 


(3-15) 


fw=*  -  xs 


Substituting  (3-15)  into  (3-14),  and  using  the  identities 


J  ifwa  —  x: 


dx  =  O 


(3-16) 


an 

f 

J  Jwa  —  x: 


dx  = 


and  CIO,  Appendix 3 


T„(x’/w)  rYk 

■  ■'■■■  ■  log)  —  Ix-x’lJ  dx 

4W*  -  x*a  l  2  J 


rrlogfYkw/4)  if  n=0 


(3-17) 


-  -  Tr,(x/w)  if  n>0 


where  T„(x)  is  the  Chebyshev  polynomial  of  the  first  kind 
and  nth  order,  (2—33)  then  becomes 


-2  log ( Ykw/4) 


r,  n 


1  -  (kw/2) 


(w*/2) (kw/2)2| |U„3 


(3-18) 


By  inspection,  then. 


12  2 

b i  =  -  -  -  log(Ykw/4)  s - log(Vkw/4) 

TT  1  -  Ckw/2)  *  n 

(3-19) 

C 

«,  =  - - 

-/w2  —  X2 


and 


8 

ba  =  — — — 
ir  (kw)  2 

(3-20) 

Dx 

"2  =  — == 

-/w2  -  X2 


are  the  solution-pairs  -for  (3-18).  In  (3-19)  and  (3-20),  C 
and  D  are  constants  to  be  determined  according  to  (2-36) . 
Hence 


(3-21) 


2^2 n/k  x 

Ma  =  -  - 


(3-22) 


Furthermore,  it  readily  -follows  from  (3-19)  and  (3-20)  that 


the  characteristic  values  of  the  narrow  slot  are  very  large 
positive  numbers  whose  ratio  satisfies 

ba  -  4  -4  4 

—  =  -  =  -  >  —  (3-23) 

b»  (kw)*  log(ykw/4)  C  (kw)  log  (ykw/4)  3  (kw)  kw 

The  equivalent  magnetic  current  M  of  the  narrow  slot 
is  now  given  by 


<M.,  I>  <Ma,  I> 

M  =  -  M»  +  -  Ma 


1  +  jb» 


1  +  jb3 


<Ma,  I>  <Ma,  I> 

- -  j -  Ma 

1  +  Jba  ba 


n 

j- 


Hy-=  (0) 


+  e«»«=(0) 


2  [  k  v+log(kw/4)  +jw/2 


*]  — 


(3-24) 


In  (3—24),  v  =  log(y)  is  Euler’s  constant.  Higher  order 
solutions  can  be  obtained  by  retaining  more  terms  in  the 
small  argument  approximations  of  Bessel  functions  in  (3-10) 
and  using  more  expansion  functions.  Incidentally,  the 
magnetic  current  (3-24)  is  identical  with  the  solution 


given  in  £103. 


3 


3-3  Evaluation  of  the  Matrices  G  and  B 

The  solution  o-f  the  eigenvalue  equation  (2-11)  -for  an 
arbitrary  slot  can  be  effected  by  transforming  the  integral 
eigenvalue  equation  into  an  algebraic  one  as  is  seen  in 
Section  2-6.  However,  the  evaluation  of  the  matrices  G  and 
B  constitutes  a  large  portion  of  the  work  involved  in  the 
solution.  An  efficient  evaluation  of  these  matrices  is 
therefore  necessary  for  the  success  of  the  solution. 

The  1 kth  elements  of  G  and  B  are  given  by 


n 

■  Gik 
k 


+w 

J 

-w 


+w 

dx  J  f  i*  (x  ’  ) 
-  w 


Jo (k I x-x 


i  )  dx’ 


n 

Blk 

k 


(3-25) 


+w 

I  <,<x’ 

-w 


+w 

dx  J  f*(x 
-w 


)  No(klx-x’l)  dx’ 


where  the  real  functions  f*  are  so  far  unspecified.  A 


particularly  simple  choice  for  f*  is 


which  is  a  pulse  expansion  -function  (see  Figure  3-2). 
Using  (3-26),  the  lkth  elements  o-f  S  and  B  become 


U  Gitt 
k 


Xx  Xk 


Jo  <k  i  x-x ' 


I ) 


X 1  —  X  X  k-1 


dx ' 


n 

—  Bi  k 

k 


Xk 

No  ( k  I  x  — x  '  I  ) 
Xk-t 


dx  ’  . 


Put 


Gjk(x) 


Xk 

Jo<klx-x,l>  dx' 

Xkt-x 


x€S*. . 


i.  (x ) 


Xk 

No  <  k  I  x -x  ’  I  )  dx  ' 

X  k —  X 


(3-27) 


(3-28) 


Then,  by  the  -first  mean  value  theorem  -for  integration  C21, 
Section  7-18 J,  there  exist  points  xJO  and  Xno(Si  such  that 


Gx  kc  =  t  Si  I  Gik  (xjo) 
k 

(3-29) 


k 


I  Sx I  Btk (Xno) . 


The  evaluation  of  GXI< 


and  61  k 


is  completed  by  integrating 


the  Bessel  functions  J0  and  N©  over  Sk,  and,  for  that 
purpose,  any  quadrature  rule  can  be  used.  Thus 

n  1  Q  p 

-  Gx*  =  -  I  Si  l  I  Si«  I  Z  q»  Jo  klxjo  -  xtk| 
k  2  i =1  L 

n  1  Q  r 

"Bit.  -  -  -  I  Si  II Sk I  Z  q±  No  klxNO  -  xt 

k  2  1=1  L 

In  (3-30),  Q  is  the  order  of  the  rule 
coefficients,  and  its  abscissas  are  given  by 

ISk  I 

—  X|<-i/2  +  Pi  — — 

2 

where  Xk-ixat  is  the  midpoint  of  Sk. 

When  evaluating  the  diagonal  elements  of  B  (l=k).  No 
offers  a  logarithmic  singularity  at  x-x’  that  requires 
particular  attention.  Put 

N0„(klx-x'l)  =  (No  -  N0» ) ( k I x -x ' I )  (3-32) 

where  No.  is  the  singular  part  of  No  given  by  its  small 
argument  approximation  (3-10).  Then 


(3-30) 


I  - 


,  q4  are  its 


(3-31 ) 


36 


*53 


a 


-  Bn  = 
k 


-I 


Xi  Xi 


dx  I  (No-  +  No— )  ( k  I  x -x '  I  )  dx’ 


*1—1  Xi-i 


,s*'’  f2Io9[Yk!!r]  -3] 


Xi  Xi 


dx  I  No- ( k I x-x ’ I )  dx ’ . 


Xi-1  Xi-1 


(3-33) 


N— „  has  no  singularity  at  x=x’,  and  can  therefore  be 
integrated  by  the  -first  mean  value  theorem  and  quadratures. 
Thus 


-BXi  = - I  Si  I*  2  log 


[Yk-rl  -3 


1  Q 

-  ISx I*  z  q4  No-  k I Xnof 

2  i  =1  L 


-  x4»|J 


(3-34) 


for  some  Xno-€Si. 


Actually,  -finding  such  points  x00,  xMO,  and  xMO-  is  at 
least  as  difficult  as  computing  the  integrals  themselves. 
For  sufficiently  small  IS*!,  however,  the  midpoint  of  Sx 
can  replace  these  points  while  introducing  negligible 
error.  Thus,  dividing  by  the  factor  ISil, 


"  G»  k 

k  I  Si  I 


-  IS*  I  £  q4  J0(Ui‘k) 
2  i  =  l 


■for  all  1 ,  k 


1  Q 

- IS*I  £  q4  Nolui11*)  if  lj^k 

2  i=l 


where 


1  Q 

-  Z  qiNop(utl,«) 

2  1  =  1 


if  1  =k 


=  k -  Xlk|.  (3-36) 

The  matrices  6  and  B  so  obtained  are  the  same  as  those 
resulting  from  enforcing  the  point  matching  condition 


R(Xi_1/a>  =  0  1=1,2, _ , N  (3-37) 

in  (2-30),  rather  than  the  Galerkin  condition  (2-31), 
except  for  a  slight  alteration  in  Bn,  as  can  easily  be 


establ i shed 


3-4  Numerical . Resul ts 

The  character!  sti  c  currents  and  magnetic  -fields  as 

well  as  the  equivalent  magnetic  current  and  radiation 

pattern  have  been  computed  -for  di-f-ferent  slot  widths.  The 
computational  aspects  of  the  solution,  and  some  of  the 
results  obtained  for  slots  of  width  0.4  A  ,  0.5A  ,  and  l.uA 
where  A  is  the  wavelength  are  discussed  in  this  section. 

In  the  actual  computation,  polynomial  appro;;  1  mati  ons 

of  the  Bessel  functions  J©  and  N©  C22,  Articles  9.4.1- 

9.4.33  are  utilized,  while  all  the  integrals  are  computed 
using  an  eight-point  Gaussian  quadrature  C22,  Table  25.43. 
In  evaluating  the  far  fields,  the  Hankel  function  is  first 
replaced  by  its  large  argument  approximation  19,  Appendix 
D3: 

38) 

C9, 

39) 

for  all  points  (x,z)fCrr,  r>>2w,  where  0  is  the  angle 
r=xu„+zu«  makes  with  the  x-axis,  the  electric  vector 


39 


potential  produced  by  the  characteristic  current  M„  at  any 
point  (r,0)€Cv<-  is  readily  -found  as 


±1 


F(M„)  = 


•<12  j  irk  r 


-jkr  N 
?  2  U„* 

k=l 


X* 


j  kx ' cos0 

e  dx  ’  u. 


N  jkx*<_t^3t  cos0 

*  ±  A(kr)  I  U„„  IS*I  e 
k=l 

sin  (  (k  IS*  Icos0)  /  2) 

-  uv  (3-40) 

(k  IS* I cos0) /2 

where 

-jkr 

e 

A ( kr )  =  - 

^2jwkr 


and  where  the  upper  sign  is  used  for  z<0  and  the  lower  one 
for  z>0.  The  far  characteri stic  fields  then  follow  by 
substituting  (3-40)  into  (3-3) . 


An  IMSL  Library  2  subroutine  "EIGZF"  C233  is  used  to 
solve  the  matrix  eigenvalue  equation  (2-33)  for  the 
characteristic  values  and  currents.  In  all  the  computer 
runs,  a  "performance  index"  has  consistently  been  less  than 
one,  indicating  that  the  subroutine  has  performed  well. 
The  convergence  patterns  for  the  characteri sti c  values  for 
the  0.4  X  ,  0.5  X  ,  and  1.0  X  slots  are  shown  in  Tables  3- 
1,  3-2,  and  3-3,  respectively.  As  can  be  seen,  the 
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The  convergence  of  the  characteristic  values  for  the  0.4  X 
slot. 


4  0.386499  2.541965  52.985104  7.204746x10  t 

3  5 

8  0.382960  2.381441  39.523594  2.778433x10  4.624070x10 

3  5 

12  0.381978  2.340359  37.074282  2.372682x10  3.329886x10 

3  5 

16  0.381576  2.322112  36.091969  2.232416x10  2.945282x10 

3  5 

20  0.381381  2.311973  35.573169  2.163472x10  2.775236x10 

3  5 

24  0.381277  2.305593  35.256428  2.123233x10  2.700021x10 
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Table  3-2 


The  convergence  of  the  characteristic  values  for  the  0.5 X 
slot. 


N 

b 

b 

b 

b 

b 

m  m  m 

1 

2 

3 

4 

5 

3 

4 

0.248621 

1 . 722857 

22. 177586 

1.848905x10 

* 

» 

3 

5 

8 

0.251433 

1.615827 

16.842234 

0.726753x10 

0.767466x10 

<D 

3 

5 

12 

0.252425 

1.588517 

15.861717 

0.622888x10 

0.554707x10 

m 

3 

5 

16 

0.252971 

1.576381 

15. 466336 

0.586800x10 

0.494704x10 

a> 

3 

5 

20 

0.253329 

1 . 569636 

15.256797 

0.569042x10 

0.466199x10 

m 

3 

5 

24 

0. 253587 

1.565691 

15. 128519 

0.558653x10 

0.450576x10 

<D 

convergence  is  always  monotone,  either  upward,  or  downward. 
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For  the  1.0  A  slot,  ba  -first  decreases  monotoni cal  1  y  and 
then  increases,  but  this  can  be  attributed  to  rounding 
errors.  Also,  the  convergence  of  the  lower  order 
character i sti c  values  is  generally  faster  than  that  of  the 
higher  order  ones. 

Any  calculated  value  of  bn  is  probably  inaccurate 
whenever  tb„/bi  I  >10,>,t*  where  NS  is  the  number  of 
significant  figures  retained  during  the  calculation. 
Whenever  bn  was  large  enough  to  be  unreliable,  it  was 
arbitrarily  set  equal  to  ®.  Now,  the  contribution  of  these 
currents  is  extremely  small,  and  therefore  has  not  been 
considered  in  subsequent  computations.  It  appears  that  for 
a  slot  in  a  conducting  plane  in  an  unbounded  medium,  only  a 
finite  number  of  characteristic  currents  need  to  be 
computed.  This  is  also  expected  to  carry  over  to  slots  in 
a  conducting  plane  separating  contrasting  mediums. 

The  computed  characteristic  currents  normalized  to  a 
maximum  amplitude  of  unity  and  their  radiation  patterns  for 
the  slots  considered  are  shown  in  Figures  3-3,  3-5,  3—7,  3- 
9,  3-11  and  3-13.  Figures  3-5,  3-9,  and  3-13  are  polar 
plots.  The  normal i zati on  to  a  maximum  amplitude  of  unity 
is  used  only  for  plotting  convenience.  The  equivalent 
magnetic  currents  and  radiation  patterns  for  the  slots  are 


shown  in  Figures  3-4,  3-6,  3-8,  3-10,  3-12,  and  3-14. 
Figures  3-6,  3-10,  and  3-14  are  polar  plots.  All  currents 
clearly  exhibit  the  right  behavior  at  the  edges.  It  is 
interesting  to  note  that  the  number  of  lobes  in  each 
pattern  is  equal  to  the  order  of  the  characteristic  current 
or  field.  When  a  slot  is  excited  by  the  plane  wave  (3-1) 
with  8=0°,  the  equivalent  magnetic  current  M  is  given  by 
(2-15).  The  power  radiated  by  the  magnetic  current  Vr,M0  is 
called  P„.  For  various  values  of  n,  the  ratio  Pn/Px  is 
given  in  Table  3—4.  This  ratio  was  evaluated  with  N  as 
specified  in  the  arrow  marked  rows  in  Tables  3-1,  3-2,  and 
3-3.  The  entries  in  Table  3-4  suggest  that  the  radiation 
pattern  for  the  slot  is  basically  the  same  as  that  for  its 
dominant  character i st i c  current.  This  is  indeed  the  case, 
as  is  readily  established  by  comparing  the  corresponding 


figures 


Table  3-4 


The  ratio  o-f  the  power  radiated  by  each  characteristic 
current  to  that  radiated  by  the  dominant  characteristic 
current  -for  a  (a)  0.4 A  slot,  (b)  0.5  A  slot,  and  (c)  1.0 A 
slot. 
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1.0 
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n 

p  /  p 
n  1 
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1.0 

2 

0.267814 

3 

-2 

1. 145190x10 
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( c ) 


n 

P  /  P 
n  1 

1 

1.0 

2 

0. 149602 

3 

0.039083 

O' 

4 

3. 766960x10 

-4 

5 

1.726436x10 

-6 

6 

5.240201x10 

-7 

7 

1.126339x10 

a  a  a 

0.  0 

radiation  patterns  for  the  first  5  i 
c  magnetic  fields  for  the  0.4X  slot. 


radiation  pattern  for  the  0.4X  slot 


e  first  5  characteristic  currents  for  the 
5A  slot. 


equivalent  magnetic  current  for  the  0.5A  slot 


magnetic  fields  for  the  0.5X  slot. 


equivalent  magnetic  current  for  the  1.0A  slot 


magne 
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CHAPTER  4 

TRANSVERSE  MAGNETIC  (TM)  CASE 

The  case  of  a  uni-form  transverse  magnetic  (to  the  slot 
axis)  plane  wave  excitation  o-f  the  slot  is  considered  in 
this  chapter.  The  character i st i c  currents  and  electric 
-fields  as  well  as  the  equivalent  magnetic  current  and 
radiation  pattern  are  determined  for  different  slots. 
Analytic  expressions  for  the  narrow  slot  are  also  given. 


4— 1  Basic  Formulation 

Let  a  plane  wave  be  incident  on  the  slot  at  an  angle  0 
from  the  left  of  the  screen  (see  Figure  4—1).  This  wave  is 
assumed  uniform  and  transverse  magnetic  to  the  y-axis,  and 
therefore  has  the  field  distribution 

-jk(x  sine  +  z  cosO) 

E*  =  e  Uy 

(4-1 ) 

1  r-  -l  -jk(x  sine  +  z  cose) 

H‘  =  -  -  | cose  u„  -  Sine  u«  e 

Since  the  slot  is  uniform  along  the  y-axis,  and  since 
the  incident  electric  field  has  only  an  Ey  component  that 
does  not  vary  with  y,  so  does  the  scattered  electric  field. 
Consequently,  the  total  field  is  independent  of  y  and 


transverse  magnetic  to  the  y-axis.  It  then  -follows  -from 
(2-1)  that  M  has  only  an  x-component  that  does  not  vary 
with  y: 


M  =  M(x)  u„ . 


(4-2) 


The  field  to  the  le-ft  o-f  the  screen  is  thus  given  by 
C9,  Section  3-123 

E(M)  =  -  V  x  F(M> 

(4-3) 


1 

H  (M)  =  -  j —  CV  (V  .  F  (M)  )  +  k2  F  (M)  3 
nk 


plus  the  short  circuit  -field 

-jkx  sin© 

§•«=  =  -  j  2  sin(kz  cosd)  e  uy 

(4-4) 

2  r 

H-es  =  -  -  jcos©  cos(kz  cos©)  uM  + 

-j  -jkx  sin© 

j  sin©  sin(kz  cos©)u,  e 

In  (4-3),  F(M)  is  the  electric  vector  potential  produced  by 
M  in  the  presence  of  the  complete  screen,  viz.  C9,  Section 
5-73: 


+w 

1  (2) 

F (M)  =  —  M(x’)  Ho  (kTz2  +  <  x  -x  ’  )  2 )  dx '  u„ .  (4-5) 

2j  , 

-w 


The  field  to  the  right  of  the  screen  is  then  the  negative 
of  that  in  (4—3) . 


Substituting  (4-3),  (4-4),  and  (4-5)  into  (2-5),  Y(M) 

and  I  are  readily  found  as 


+w 


1  d2 

Y  (M)  =  —  ( -  +  k2) 

k  n  dx2 


(2) 

M(x’)  Ho  (klx-x’l)  dx 


-w 


(4-6) 


-jkx  sin6 


I  =  -  -  cos©  e 


where,  as  before,  scalar  quantities  are  used  since  Y(h)  and 
I  are  both  directed  along  the  x-axis.  Furthermore,  using 
(3-7),  the  operators  G(M)  and  B(M)  are  now  given  by 


+w 


1  d2 

G(M)  =  —  ( -  +  k2) 

k  H  dx2 


M  (x ’ )  Jo (k I x— x '  I )  dx 


-w 


(4-7) 


+w 


1  d2 

B  (M)  = - ( -  +  k2) 

kH  dx 2 


M(x’)  N0(klx-x!'l)  dx 


-w 


Below,  the  narrow  slot  is  treated  first,  and  is  then 
followed  by  the  general  case  of  an  arbitrary  slot. 


4-2  The  Narrow Slot 

Approximations  similar  to  those  utilized  in  the  TE 
case  can  be  used  to  simpli-fy  the  analysis  of  a  narrow  slot 
excited  by  a  uniform  transverse  magnetic  plane  wave. 
Specifically,  the  Bessel  functions  Jo  and  No  are  replaced 
by  their  small  argument  approximations  (3-10). 

Furthermore,  for  a  narrow  slot,  I  does  not  vary  greatly 
over  the  slot,  and  can  therefore  be  approximated  by  the 
first  two  terms  of  its  Taylor  expansion  about  the  origin, 
viz.: 

I  ~  H-*=  (O)  +•  x  H-=*  (0) 

=  -  -  ^2cos8  -  jkx  sin (20) j.  (4-8) 

In  (4—8),  the  prime  denotes  differentiation  with  respect  to 
x,  and  (0)  is  written  for  the  point  (x,y)=M0,0>. 

The  characteristic  currents  are  determined  by  solving 
the  matrix  eigenvalue  equation  (2-33).  Since  the  slot  is 
narrow,  however,  partitioning  of  it  is  not  necessary.  The 
lkth  elements  of  G  and  B  then  become 


Gx »« 


k y\ 


j 


f  i  ( :<  ) 


-w 


dx 


d= 

dx2 


+w 

2 

1 

r  k  -i 

fh  (X'  ) 

1  -  — lx-x ' I2 

l 

L  4  J 

dxJ 


+w 


B 1  •«  — 


-2 

=  -  f  i  (x 

wkn  , 


)  dx 


-w 


(4-9) 


+w 

d*  P  rYk 

( -  +  k2) I  f„(x'>  log!  —  Ix-x 

dx2  J  1-2 

-w 


An  appropriate  choice  of  f  that  satisfies  the  edge 
requirements  indicated  in  Section  2-2,  and  is  compatible 
wi th  (4—8) ,  is 


f * ( x )  =  x*~l  ^w2  -  x2  k=l , 2 


(4-10) 


which  results  in  a  second  order  matrix  eigenvalue  equation. 
Substituting  (4-10)  into  (4-9),  and  using  the  identities 


n  >0 


+w 


(2n-l )  ! 


w 


2 (n^l > . 


x2”  -  xa  dx  = 


(n+1 )  !  (n-1 )  !  22" 


-w 


and  CIO,  Appendix] 


+w 

at 

f  rYk 

i  w  IT  rYkw-1 

|  -/w2  -  x  ' 2  log! —  1  x— x  ’  1 

dx’  =  -  log  - 

J  L  2 

J  2  L  4  J 

— w 


w  n  r 

+  -  |2(x/w> 


-*] 


(4-12) 


+w 


j*  Ur,(x’/w)  -Jw2  -  x '  *  log  j-—  I  x  -x  ’  I  j  dx’  = 


-w 


w2tr 


n (n+2) 


JVl  J^l  -  <X/w)2j  Ur, 


(x/w)  +  Tr»  (  X  /  W  ) 


]  n>0 


where  T. 


and  Un 


respectively. 


the  Chebyshev 


are, 

polynomials  of  the  first  and  second  kinds,  <2-33>  then 
becomes 


r  rYkw-. 

r 

2  -  (kw/2)2  +  (kw)2log -  0 

L  4  J 

r  2  1 

0  w2  1-  -(kw/2)2 

U„a 

L  L  3  JJ 

. 

r  r  i 

■ 

•  ' 

(kw/2) 2 1 1  -  (kw/2) 2  I 

0 

U„t 

w2ir2 

L  J 

b  - 

"  2 

0 

(w/2) 2 (kw/2) * 

Ui2 

. 

■  . 

n=l ,  2 
(4-13) 


By  inspection 


-1 


b,  = 


ir(kw/2)2Cl  -  (kw/2) *3 


2  -  (kw/2) *  +  (kw)2log 


b 


(4-14) 


M»  (x)  *  C  -Jw2  -  xs 


and 
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are  the  solution-pairs  of  (4-13).  As  can  be  seen,  b»  and 
b2  are  very  large  negative  numbers.  Furthermore,  since 
2kw<<l,  the  characteri sti c  values  are  well  approximated  by 


4  j-Y  kM-< 

-  -  2  +  ( kw)  21  og  I - | 

tt  (kw)  a  l  4  J 


(4-16) 


it  (kw) 


Consequently,  to  this  approximation,  the  characteristic 


(4-17) 


currents  are  given  by 


2-TZtt 

M»  (x)  =  -  -Jw3  -  x! 

kW=*7T 


8-/2r)k 

Ma(x)  =  -  X  -Jwa  -  Xs 

kaW^TT 


when  normalized  according  to  (2-36). 


The  equivalent  magnetic  current  on  the  slot  is  given 


<Mi,  I>  <Ma,  I> 

li  (x )  =  -  Mi  +  -  Ma 

1  +  jb»  1+jba 
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< M x  ,  I>  <M»,  I>. 

-  M*  -  j -  Ma 

1  +  jbx  b* 


w  TJ  k 


H-«=(0> 


(kw/2)2  —  j  <2  +  (kw)alog  (y  kw/4)  ) 


-  H-«*  (O)  xj  iw2  ->:2. 


(4-18) 


The  magnetic  current  (4-18)  reduces  to  that  in  CIO!  when 
all  terms  containing  kw  -factors  are  ignored.  Higher  order 
solutions  can  be  obtained  by  retaining  more  terms  in  (3-10) 
and  using  more  expansion  -functions. 


4-3  Evaluation  of  the  Matrices  G  and  B 

The  general  slot  of  arbitrary  width  can  be  solved 
using  the  Galerkin  procedure  of  Section  2-6.  Similar  to 
the  TE  case,  an  efficient  evaluation  of  the  matrices  G  and 
B  is  essential  for  the  success  of  the  solution. 

Because  of  the  form  of  the  operators  G  and  B,  a 
triangular  expansion  of  the  characteri sti c  currents  (see 
Figure  4-2)  is  attempted.  This  has  the  advantage  of 
reducing  the  second  order  partial  derivatives  to  finite 
differences  C24D,  as  well  as  approaching  zero  as  x-*±w, 
although  not  at  the  rate  specified  by  the  edge  conditions. 


To  keep  the  notation  simple,  a  uniform  partition  of  the 
slot  is  assumed.  Thus 


where  S  is  the  uniform  width  of  the  partition.  The  lkth 
elements  of  G  and  B  are  then  given  by 


Xi  +  t 

lk  =  f  X  ( X  )  I 

I  ** 


x  -  au-± 


-  x 


Xk-l  <  X  <  x, 


X*  <  X  <  xte*, 


( -  +  k*)Gte (X)  dx 

dx*  J 


X  K  =  -  f  1  <X  )  I 


( -  +  k*)B* (x) |  dx 

dx*  J 


where 


G„  <x) 


fk(x')  Jo(klx-x’l)  dx’ 


f*  <x’ >  N©(klx-x’t)  dx’ 


(4-19) 


(4-20) 


(4-21) 


w«5 


Integrating  the  -first  term  in  (4-20)  by  parts  twice, 
and  approximating  G*(x)  by  G* (x» )  and  B„(x)  by  B* (Kx )  in 
the  second  term,  Gx*«  and  Bx*  then  become 


>x*«  —  —  |gh(xi*i)  —  2Gi,(xi)  +  Gk(xi_t)| 


Xi*i 


f  n 


+  ka  Gi«(Xi)  I  -f  x  (x )  dx 
J 

X  X  —  X 


_ _ i  r, 

s  L 


Bi«  (Xx*t)  —  2B„  (x i  )  +  Bk  (x 


(4-22) 


-  k*  Bh  (x 


AX* 

*’  lu 


(x)  dx, 


Xx-1 


Since 


Xx*x 


•fx  (x)  dx  = 


(4-23) 


Xx-x 


the  elements  Gi*  and  BXt<  o-f  the  matrices  G  and  B  become 


knGx  i 


—  -  [ft.  (Xi*x)  — 

s  L 


’  ( 1  -  DGk(xi)  +  Gkt  (x 


->] 


(4-24) 


knBx, 


=  -  i  r 


Bt«  (Xx*x)  -  2(1  -  DB 


h(xx)  +  Bi«  ( x  i  —  t  >1 


Xi»i,  and,  -for  that  purpose,  any  quadrature  rule  can  be 


used. 


When  evaluating  the  diagonal  and  the  -first  and  second 
lower  and  upper  diagonal  elements  o-f  B,  Bik,  ll-kl<2. 


singularities  are  encountered.  As  has  been  seen  in  the  TE 
case,  this  situation  is  best  dealt  with  by  integrating  the 
singular  part  o-f  the  integrand  analytically.  The  non¬ 
singular  remainder  is  then  integrated  numerically  using  a 
quadrature  -formula.  Thus 


Bk(xj)  *  B*<xj)  +  Bh(Xj) 


P 

B*  (x 


X  K—  1 


■f*<x’)  No«»  <k  lxj-x  ”  I  )  dx 


X*-l 


. I 


X|«*t 


■flt(x’)  Noa  (  k  I  X  J  —x  '  I )  dx 


(4-25) 


where  No.  is  the  small  argument  approximation  o-f  No  given 


integration 


When 

s 

(x 


+  I  j-i  I 


rYk 

logj^ — J  j-i  IS 


(4-28) 


s 

Ij— kl=l,  Bk(xj)  o-f  (4—28)  reduces  to  (1*  and  when  j=k, 
)  reduces  to  (5*  where 


(4-29) 


The  quantities  (B*(xj>,  j=l-l,  1,  1+1)  appear  on  the 


right-hand  side  o-f  the 


second  o-f  equations  (4-24). 


Retaining  (4-24)  whenever  I  k-1  I  >2,  and  replacing  B*  (xj  )  by 
means  of  (4-25)  and  (4-29)  whenever  lj-kl<l,  we  obtain 


G*  (Xx  -*-x  )  —  2(1  —  D  Gx«  (xx  )  +  Ghc  (xi-i  ) 

«.  - 


■for  all  1  ,  k 


1  a  I- 

-  B*(> 

1  S  L 


(xl  +  l)  -  2(1  -  r)Bh(Xx)  +  B^xx-x) 


if  1 k-1 I >2 


1  r p 

-  -  Bk(x,*i)  -  2(1  -  r)Bk(Xi)  +  Bk(: 
S  L 


x-x)  +  ftxj 


if  k-l=; 


B|«  (Xx*x)  —  2(1  —  DBk(Xx)  +  Bk  (x  x  -  x )  +  ft: 


if  i-k=: 


1  r p  p 

!-  —  Bk  (Xx*i)  —  2(1  —  n  B|*  (Xx  )  +  Bk  (Xx-x  ) 
;  S  L 


-  -  k(Xi. 
s  L 


+  Ga  -  2(1  -  n  fix 


p  p 

)  -  2(1  -  r>B*(xx)  +  Bk  (xx  _x ) 


+  fta  -  2(1  -  D 


if  k— 1=1 


if  1 -k=l 


1  r  P  P  P 

—  I (Xx*i )  —  2(1  —  r)  Bit  ( x i  )  +  Bit  ( x x  —  i ) 
S  L 


+  2ft i  -  2(1  -  D  ft: 


if  1 =k . 


(4-30) 


In  (4-30),  B„  (>{j),  j€£k-l,  k,  k+13,  is  evaluated  as 


P 

B„  (Xj  ) 


S  E  c„  f  *  (xte+dr,S)  Nop  ( k  I  x  j  —  (X|,+dnS )  I  ) 
n=  1 


(4-31) 


where  Q  is  the  order,  c„  are  the  coefficients,  and  d„  are 
the  abscissas  of  the  quadrature  rule. 


4-4  Numerical  Results 


The  character i sti c  currents  and  electric  fields  as 
well  as  the  equivalent  magnetic  current  and  radiation 
pattern  have  been  computed  for  different  slot  widths.  Some 
of  the  results  obtained  for  the  0.4A,  0.5  A,  and  1.0  A  slots 
are  given  in  this  section. 


In  the  actual  computation,  polynomial  approximations 
of  the  Bessel  functions  J0  and  No  C22,  Articles  9.4.1- 
9.4.33  are  utilized,  while  all  the  integrals  are  computed 
using  an  eight-point  Gaussian  quadrature  C22,  Table  25.43. 
The  char acteristic  electric  fields  E„y  are  given  by 


±1  kz 

-nv  =  ” "  ~ 

j 2  -Jza+(x-x’  ): 


Mr,  ( x  ’  )  Ht  (kJz2+(x-x’  )2  )  dx’ 


(4-32) 


# 


is  the  Hankel  -function  of  second  kind  and  first 


where  Ht 

order.  The  upper  sign  is  used  for  z<o  and  the  lower  one 
for  z>0.  The  far  fields  can  then  be  obtained  by  replacing 
the  Hankel  function  by  its  large  argument  approximation  [9, 
Appendix  DI: 

(2)  _  -js 

Hi  <  5 )  =  j  -Jf2j  /  (tts)  e  (4-33) 

and  using  the  radiation  zone  approximations  (3-39).  Thus, 
substituting  (2-28)  for  M„  with  f„(x)  given  by  (4-19),  we 
obtai n 


1+j 


E«w  X 

ny  ~  — 


Jk7  (irr)  sin0 


|-sin  (  (kScos0) /2)-j2  -jkr 

- j  s  e 

l  (kScos0)/2  J 


L  jkxkcos0 

E  U„k  e 
k=l 


L  jkxt<cas0 

*  A (r , 0)  Z  U„*  e  (4-34) 

k=l 

at  any  point  (r,0)€Cy,r.  when  r>> 2w  and  kr>>l.  In  (4-34),  0 
is  the  angle  r=xuK+zu,  makes  with  the  x-axis. 

All  the  features  and  convergence  patterns  found  for 
the  TE  case  are  repeated  for  the  TM  case.  The  convergence 
of  the  character i sti c  values  is  monotone,  either  upward  or 
downward,  as  can  readily  be  seen  in  Tables  4-1,  4-2,  and  4- 


The  convergence  of  the  characteri sti c  values  for  the  o.4A 
slot. 


Table  4-2 


The  convergence  o-f  the  characteristic  values  -for  the  o. 
slot. 


L 

b 

b 

b 

b 

1 

2 

3 

4 

4 

-0.4591318 

-14.9436836 

-897. 567627 

<D 

8 

-0.4725091 

-14.0017996 

-648.799561 

-62683.6680 

12 

-0.4709589 

-13.6431913 

-600.401367 

-51791.5703 

16 

-0.4687450 

-13.4430857 

-579.949951 

-49074. 1836 

20 

-0.4668146 

-13.3114119 

-567.934814 

-47766.3203 

24 

-0.4653491 

-13.2199717 

-560.609375 

-48250.5898 

28 

-0.4641726 

-13. 1521626 

-555.927002 

-49039. 1602 

Table  4-3 


The  convergence  o-f  the  characteristic  values  -for  the  1.0  A 
slot. 


4  0.007760  -0.2222299  -6.3792248  -316.052002 


8  -0.0100392  -0.4589630  -7.8628264  -187.909714  -8264.32031 


12  -0.0139949  -0.5044749  -8.1955252  -179.558655  -6616.984337 


16  -0.0157180  -0.5185139  -8.2802639  -176.418991  -6160.37109 


20  -0.0167178  -0.5237506  -8.2951450  -174.405258  -5933.50781 


24  -0.0173952  -0.5258325  -8.2879648  -172.971313  -5828.77734 


28  -0.0178963  -0.5265995  -8.2730077  -171.782654  -5732.35547 


32  -0.0182855  -0.5267935  -8.2571888  -170.876831.  -5652.25781 


36  -0.0185869  -0.5265971  -8.2394047  -170.050171  -5601.01953 


40  -0.0188513  -0.5263984  -8.2251053  -169.446167  -5562.56641 


3.  It  should  be  noted,  however,  that  only  four 
characterl st i c  modes  for  the  O.AX  and  0.5  A.  slots  and  five 
for  the  1.0  A  slot  are  computed  compared,  respectively,  to 
five  and  seven  for  the  TE  case.  This  is  because  the 
corresponding  programs  were  run  on  different  computers 
having  different  computing  precisions.  When  the  TE  program 
was  run  in  the  new  computing  environment,  there  resulted 
the  same  number  of  significant  modes.  The  characteristic 
currents  and  far  electric  fields  for  the  slots  considered 
are  shown  in  Figures  4-3,  4-5,  4-7,  4-9,  4-11,  and  4-13. 
Figures  4-5,  4-9,  and  4-13  are  polar  plots.  The  curve  of 
the  characteristic  current  l"U(x)  in  Figure  4-3  is  not 
exactly  odd  about  x=0  because  of  roundoff  error.  The 
corresponding  equivalent  magnetic  currents  and  radiation 
patterns  are  shown  in  Figures  4-4,  4-6,  4-8,  4-10,  4-12, 
and  4-14,  Figures  4-6,  4-10,  and  4-14  are  polar  plots. 
When  a  slot  is  excited  by  the  plane  wave  (4-1)  with  8=0°, 
the  equivalent  magnetic  current  M  is  given  by  (2—15). 
Table  4-4  gives  for  each  slot  the  ratio  of  the  power 
radiated  by  each  characteristic  current  in  (2-15)  to 
that  radiated  by  the  dominant  characteristic  current  V*M». 
All  data  are  evaluated  with  L  as  specified  in  the  arrow 
marked  rows  in  Tables  4-1,  4-2,  and  4-3. 
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Table  4-4 

The  ratio  o-f  the-  power  radiated  by  each  characteristic 
current  to  that  radiated  by  the  dominant  characteristic 
current  -for  a  (a)  0.4  X  slot,  (b)  0.5  X  slot,  and  (c)  1.0  X 
slot. 


m 


n 

P  /P 

n  1 

1 

1.0 

2 

0. 100758 

-2 

3 

0.251782x10 

-4 

4 

0. 1633407x10 

•  •  • 

O 

« 

o 

(b) 

n 

p  /p 

n  1 

1 

1.0 

2 

0. 161499 

-2 

3 

0.6648028x10 

-4 

4 

0.7618671x10 

— 

O 

• 

o 

,•  1 1  *•  , 


equivalent  magnetic  current 


ivalent  magnetic  current  for  the  1.0X  slot 


CHAPTER  5 


DISCUSSION 


The  theory  of  characteristic  modes  for  slots  has  been 


applied  to  the  solution  of  the  problem  of  an  infinitely  long 


slot  in  a  ground  plane  illuminated  by  either  a  transverse 


electric  or  a  transverse  magnetic  plane  wave.  The  theory  is 


a  specialization  of  the  general  theory  of  characteristic 


modes  for  apertures,  and  can  be  applied  to  slots  in  a 


conducting  plane  separating  contrasting  mediums  and 


illuminated  by  a  general  wave. 


Some  conclusions  are  drawn  from  the  application  of  the 


theory  in  the  report. 


1-  The  theory  is  equally  applicable  to  narrow  as  well  as 


to  general  slots. 


2-  The  theory  results  in  a  modal  expansion  for  all 


quantities  and  parameters  of  interest  usually  encountered  in 


electromagnetic  compatibility  problems,  such  as  the 


equivalent  magnetic  current  on  the  slot,  the  transmission 


coefficient,  and  the  radiation  pattern. 


3-  For  narrow  slots,  the  theory  reduces  to  an  augmented 


Bethe  theory,  i.e.,  the  slot  is  represented  by  a  susceptance 


term  related  to  the  polarizability,  plus  a  conductance  term. 


\t  v  ',1  v  y.ylty. 


«■.  f,  »»,  r 


4-  Only  a  finite  number  of  characteristic  currents  and 
field  modes  need  to  be  computed  for  a  slot  of  arbitrary 


width. 


5—  When  the  slot  is  narrow,  only  a  few  terms  are  needed 
in  the  modal  expansion  mentioned  in  2  above. 


Although  the  presentation  is  confined  to  a  slot  in  a 
conducting  plane  in  an  unbounded  medium,  the  conclusions 
drawn  above  are  expected  to  hold  for  a  slot  in  a  conducting 
plane  separating  mediums  with  different  electromagnetic 
properties.  The  extension  to  the  two-medium  problem  can  be 
accomplished  in  a  straightf orward  manner.  The  application  of 
the  theory  to  another  kind  of  aperture  like  the  circular  or 
elliptic  aperture,  or  for  that  matter  a  general  aperture  of 
arbitrary  shape,  however,  is  a  major  task  worthy  of 


consideration. 


The  study  of  the  transmission  of 


electromagnetic  energy  through  apertures  backed  by  resonant 
cavities,  or  loaded  by  dielectric  sheaths  or  strips  is  of 
both  theoretical  and  practical  importance.  Other  areas  of 
possible  study  include  phased  array  antennas,  pattern 
synthesis,  and  supergain  array  design. 


The  theory  of  characteristic  modes  for  apertures  is 
versatile.  It  has  been  applied  with  great  success  to  slots 


VMV. 


Hf  ■ 


v  ..‘N'Jv’VyVV 
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Appendix  A 

In  this  Appendix,  it  is  shown  that  the  convergence  of 
the  modal  expansion  o-f  E(H)  is  in  a  least  squares  sense  over 
the  radiation  cylinder.  The  proof  for  H(M)  is  similar.  Let 

(N)  N 

E  <M>  =  E  V„  E„  ( A— 1 ) 

n=  1 

be  the  N-term  modal  expansion  of  E(M>  and  let  E(M)  be  the 
quantity  toward  which  it  converges. 

N 

E(M)  =  lim  ZVr.gr,  (A— 2) 

N-k»  n=l 


If 


U„  -  lim  V„ 
N-»® 


(A— 3) 


then 


E(M)  =  Z  U„  E„.  (A— 4) 

n=l 


The  difference  between  <A-1>  and  (A-4)  is  called  the 
residual  R^. 


Rn  —  Z  Ur.  En 
n=  1 


N 

Z  V„  E„ 
n=  1 


(A— 5) 
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Then,  the  convergence  of  E(M)  in  a  least  squares  sense  over 
the  radiation  cylinder  is  equivalent  to  the  requirement  that 
the  square  o-f  the  norm  o-f  the  residual 


I  I Rn I i  = 


* 

—  Rim  « 


Rn  dr 


(A— 6) 


be  minimum.  Substituting  (A-5)  into  (A-6)  and  using  the 
-first  orthogonality  relationship  of  (2-25),  there  results 


H  Rn  H  =  n  Z  I  Ur 
n=l 


-  Vr,«  +  Z  ILU  t 
n=N+l 


(A— 7) 


The  quantity  (A-7)  is  minimum  when 


V„  =  U„  ,  n=l , 2, .  .  .  ,N 


(A— 8) 


As  given  by  (2—19),  V„  does  not  depend  on  N.  It  is  now 
evident  from  (A— 3)  that  V„  is  given  by  (A-8) .  Therefore  the 


modal  coefficients  V„  minimize  (A-7)  so  that  the  convergence 
of  the  modal  expansion  of  E(M)  is  in  a  least  squares  sense 


over  the  radiation  cylinder. 
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